Formation of clean dimers during gas-source growth of Si(OOl) 
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Elevated temperature STM measurements have shown that one key phase during gas-source ho- 
moepitaxy of Si(001) is the formation of clean Si ad-dimers from hydrogenated ad-dimers, though 
the mechanism for this formation is unknown. We present ab initio density functional calculations 
designed to explore this mechanism. The calculations show that there is a pathway consistent with 
the experimentally observed reaction rates, which proceeds via a meta-stable intermediate, and is 
effectively irreversible. This result fills a vital gap in our understanding of the atomic-scale details 
of gas-source growth of Si(001). 

PACS numbers: 68.43.Bc;81.10.Aj;81.15.Aa;31.15.Ew 



I. INTRODUCTION 

Gas-source growth of Si(OOl) using hydrogen-based 
precursors (such as SiH4, silane, and Si2Hg,p(iisilane) is 
of great scientific and technological interestBHfla — in 
particular, hydrogen can act as an effective surfactant, 
and has been shown to reduce roughness and intermixing 
during growth of Ge/Si alloys and pure Ge on Si(001)u. 
Understanding the reactions that occur and the interme- 
diate structures that are formed during this growth will 
enable greater control of surface and interfaces during 
growth. STM observations of the growth of Si(001) frc 
disilanc, both at room temperature following anneals 
and at elevated temperatureffiLalong with careful elec- 
tronic structure calculationsouu have mapped out the 
growth pathway. A key observation in this pathway is 
that the islands formed during growth are clean, while 
the substrate remains covered with a certain amount of 
hydrogentj. The fundamental building block in gas-source 
growth is the clean ad-dimer (as opposed to solid-source 
growth, where fast-moving ad-atoms are keyQ); yet, the 
mechanism to form such clean dimers from the hydro- 
genated dimers that occur naturally during gas-source 
growth is unknown. In particular, they are observed 
to form at 450K while desorption from the monohy- 
dride phase occurs at 790K, indicating that their for- 
mation must be completely different to the desorption of 
hydrogen from the monohydride phase. In this paper, 
we present a first-principles investigation of the mecha- 
nism for formation of clean ad-dimers from hydrogenated 
dimers, with the aim of explaining how these form at a 
comparatively low temperature. 

Disilanc (which is used in preference to silane as it de- 
composes more easily) adsorbs on Si(OQl) as SHI3 (which 
soon breaks down to form S1H2) or SiH^u, sometimes with 
accompanying hydmgen. These SiH2 groupsQ start to 
diffuse at 400-500KH When two groups are on adjacent 
dimer rows, they react to form a hydrogenated ad-dimcr 
(that is, an ad-dimer with both dangling bonds saturated 
with hydrogen, illustrated in Fig. |l] (a)) over the trench 
between the dimer rowsfa. This then decomposes to form 
clean|-ad-idimers and hydrogen on the surface at around 
450KElBu, via a pathway to be investigated in this paper. 



A hydrogenated ad-dimer (which is the starting point) 
is illustrated in Fig. [l] (a), along with a partially hydro- 
genated ad-dimer (the result of the first part of the path- 
way) in Fig. [l] (b) and a clean ad-dimer (the final point} 
in Fig. [I] (c). Once formed, the clean ad-dimers diffuseQ 
and form a square feature, which is believed to be the 
precursor to dimer stringsQ, followed by short strings of 
dimerso which later increase to form larger islandsBu. 

The calculations to be presented are based on den- 
sity functional theory (DFT) in the generalized gradi- 
ent approximation (GGA), with a plane wave basis set 
and pseudopotentials. We have searched for possible 
pathways both by applying constraints to specific atoms 
(for instance constraining a hydrogen to lie in a given 
plane) raed by using the nudged elastic band technique 
(NEB)li30 which allows accurate determination of reac- 
tion barriers given an initial approximation to a pathway. 
One key result is that the dehydrogenation proceeds via 
a meta-stable intermediate state (this is discussed fully 
in Section III and illustrated in Fig. ||). 

The rest of the paper is organised as follows: the 
next section gives details of the computational techniques 
used; this is followed by a detailed discussion of the struc- 
ture of the meta-stable state which plays a key role in the 
dehydrogenation; the diffusion pathways are then pre- 
sented, looking at the mechanism for both hydrogens, 
followed by a conclusion section. 



II. COMPUTATIONAL DETAILS 



The theory underlying DFTliallj and their application 
to electronic structure calculations have been extensively 
reviewedt-3, as has the use of pseudopotential and plane- 
wave techniquealj. The calculations in this paper were 
performed using the VAS£Lcodeli3, using the standard 
ultra-soft pseudopotentialsEJ that form part of the code. 
The approximation we use for exchange-correlation en- 
ergy is the generalised-gradient approximation (GGA) 
due to Perdew and Wang (PW91)LL3'ti We chose the 
GGA rather than the local density approximation (LDA) 
rather deliberately. As the barriers that we will be calcu- 
lating are sensitive to bonding and stretched bonds, and 



FIG. 1: Structures of: (a) the starting point, with a hydrogenated ad-dimer; (b) the end point for the first diffusion event with 
one hydrogen on the substrate; and (c) the final point with a clean ad-dimer and both hydrogens on the substrate. 



the GGA is known to be rather more accurate in these 
situations (LDA generally overbinds), we considered its 
use to be essential for this work. 

We use periodic boundary conditions, as is standard 
for DFT calculations with plane-waves, and we there- 
fore used a periodic slab for the surface, with a vacuum 
layer between the slabs. Our simulations were performed 
within a unit cell two dimers long and two dimer rows 
wide, with five layers of Si (the bottom of which was ter- 
minated in hydrogen and constrained to remain fixed in 
bulk- like positions). The vacuum gap of 6.9 A is equiva- 
lent to five atomic layers of Si, and provides sufficient iso- 
lation between vertical periodic images. We used a plane 
wave cutoff of 150 eV and a 2 x 2 x 1 Monkhorst-Pack k- 
point mesh. All these parameters were tested, and found 
to converge energy differences to better than 0.01 eV. 
The system contains an even number of electrons, but 
has various saddle points which might involve unpaired 
electrons, so we checked the effect of performing spin- 
polarised calculations for these points. The effect was 
found to be negligible (both on energies and geometries) 
and so was not used in the calculations. 

To investigate the diffusion pathways, we used two 
techniques: first, constraining the diffusing hydrogen to 
lie in a particular plane, and calculating static ener- 
gies for different locations of ,tb/? .hydrogen; second, the 
Nudged Elastic Band methoclHllil. This second method 
requires the simultaneous relaxation of a number of im- 
ages of the system, which can be done in parallel. How- 
ever, this has the potential to become extremely com- 
putationally intensive, which is why we chose to use the 
smallest realistic unit cell (with four dimers in the sur- 
face) . 

The initial exploration of the system used the first 
method (static calculations, constraining the hydrogen). 
It was using this method that the meta-stable state (dis- 
cussed in Sec. Ill) was found, and it is unlikely that it 
would have been found using the NEB without signifi- 
cant effort (for instance performing simulated annealing 
on the initial images of the system) or using ajpore com- 
plicated technique such as the dimer methodEl. The dif- 
fusion barriers presented in the paper were all calculated 
using a variant of the original technique which Actively 
seeks the saddle point, the climbing image NEBcil, with 
8 images relaxed in the chain. 

While we have calculated the diffusion barriers, we 




FIG. 2: The structure of the meta-stable state which provides 
a low energy pathway for dehydrogenation. The diffusing hy- 
drogen is bonded to a substrate dimer (A), which has broken 
one bond to a second-layer silicon (B). The ad-dimer (C) is 
now partly clean, and has formed a bond to the second-layer 
silicon (B) left by the substrate dimer. 



have not calculated attempt frequencies, which have been 
assumed to be 10 13 sec -1 , typical for diffusion processes. 
While DFT-GGA is sufficiently accurate to calculate re- 
action barriers to within 0.1 eSL-.it is not able to predict 
attempt frequencies accuratelyE3; for instance, in previ- 
ous work on solid-source growth of Si(001), it was shown 
that a factor of five in the attempt frequency was re- 
quired tOj-tuiderstand the results, but was not accurately 
predictedc3. 



III. THE META-STABLE INTERMEDIATE 

The lowest energy diffusion pathway, and the only one 
which has an energy barrier which is in line with the 
temperature at which the dehydrogenation is observed 
to occur, proceeds via a meta-stable intermediate. This 
is an unusual and rather important structure, and will 
be discussed in detail in this section. It is illustrated in 
Fig. |[ and should be contrasted with the hydrogenated 
ad-dimer illustrated in Fig. [j] (a). 

The atom labelled 'A' in Fig. | is one of the substrate 
dimers, to which the hydrogenated (and clean) ad-dimer 
is bonded; 'B' is a second layer atom in the substrate, 
to which the substrate dimer is normally bonded (for 
instance in Fig. [I] (a)); and 'C is the ad-dimer atom 
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itself, which is now clean (having started hydrogenated). 
The hydrogen is now bonded to 'A', which has broken 
its bond to 'B', while 'C has formed a bond to 'B' (not 
easily seen, owing to the geometry). 

In terms of the bonding of the atoms, the atoms A, 
B and C are all saturated as they are in Fig. [l] (a) - 
the bonding has merely cycled around (so that the A-B 
bond is now an A-H bond and a B-C bond, while the 
C-H bond is now a C-B bond) . It is this saturation that 
gives the structure its stability. While some of the bond 
angles are rather strained (in particular, the bonds asso- 
ciated with B and C make 60° angles) the bond lengths 
are all close to equilibrium, and there are no further bro- 
ken bonds, leading to an energy difference of 0.57 eV rel- 
ative to the starting point, but no more. It is interesting 
to note that there are other structures where 60° bond 
angles are found during growth of Si(001), which _also 
exhibit stability which might seem counter-intuitiveo. 

In terms of the formation of this structure, as we shall 
see in the next section, there is not a large barrier. The 
H never has to move a long distance from either A or 
C, leading to relatively strong bonds being present at all 
times; the second-layer atom B moves up slightly; and 
while the substrate dimer atom A and the ad-dimer atom 
C do move up and down respectively, they do this gradu- 
ally while maintaining their bonding. It is this relatively 
small perturbation on the overall structure, and the ease 
with which it is reached, which allows the formation of 
this state, and gives it its importance. 



IV. DIFFUSION PATHWAYS 

In this section, we describe the diffusion pathways that 
we have explored with DFT calculations. For simplicity, 
and because it is likely to be physically realistic, we al- 
low the hydrogens on the ends of the dimer to diffuse off 
independently — i.e. we consider the diffusion off one 
end of the ad-dimer while the other hydrogen remains on 
the ad-dimer. Then we allow the remaining hydrogen to 
diffuse off the now partially-hydrogenated ad-dimer onto 
the substrate. In order to avoid the complications of spin 
and half-filled bands, we maintain both hydrogens in the 
unit cell at all times (the first hydrogen to diffuse off stays 
on the substrate, illustrated in Fig. 0(b)). The three sta- 
ble points of the process (fully hydrogenated ad-dimer, 
partially hydrogenated ad-dimer with a hydrogen on the 
substrate and clean ad-dimer with both hydrogens on the 
substrate) are illustrated in Fig. |l|. The atomic positions 
during the diffusion pathways are presented below in an 
aggregated form (due to space constraints): only the po- 
sition of atoms which move significantly are shown. All 
of the atomic structures at each .step in all the diffusion 
pathways are available elsewhereEil. 

The experimental data that we are comparing against 
comes from two separate experiments: first, where the 
Si(OOI) surface was exposed to a dose of disilane (Si2H 6 ), 
annealed at different temperatures for different times, 
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FIG. 3: A graph showing energy barriers for the diffusion for 
the first hydrogen off the ad-dimer. The energies are given rel- 
ative to the starting point, while the x-axis gives the distance 
from the starting point of the hydrogen. The open circles 
show direct diffusion (proceeding without the meta-stable in- 
termediate state). The open squares show the diffusion via 
the meta-stable state (whose position is marked with a ver- 
tical dotted line at 3. 35 A, and whose structure is shown in 
Fig. [^). The lines (long dashes for direct diffusion and dash- 
dotted for diffusion via the meta-stable state) are spline fits 
to the data, and are given as guides to the eye. 



and then observed at room temperature in STlVliiita; sec- 
ond, where an elevated-temperature STM was used to 
observe the results of dosing with disilane in real time 
at different temperaturesou. The results of both these 
types of experiment are identical: around 450K, clean, 
non-rotated dimers are formed over the trench between 
dimer rows. In other words, the monohydride dimers lose 
their hydrogen to the substrate in a matter of minutes at 
this temperature (for instance, an anneal to 470K for two 
minutes led to the dehydrogenation of all ad-dimersa). 
Assuming an attempt frequency of I0 13 sec -1 and a suc- 
cessful dehydrogenation rate of 1/60 Hz, we obtain a bar- 
rier of 1.28 eV. This changes by about 0.03 eV if the rate 
is doubled or halved, giving us a good estimate of the 
likely reaction barrier. 



A. The First Hydrogen 

There are two diffusion paths considered for the first 
hydrogen diffusing off the ad-dimer: a direct diffusion 
path; and diffusion via the meta-stable considered in 
Sec. Ill and shown in Fig. |[ We will discuss these sep- 
arately, starting with the direct diffusion, and then con- 
trast their results. 

The diffusion barrier for direct diffusion is shown in 
Figure^, with open circles and long dashes. The barrier 
is 1.93 eV, which is extremely high; the reason for this can 
be seen from the atomic positions, which are illustrated 
in Figure [J. At the saddle point, the ad-dimer bond is ex- 
tended greatly (from 2.51 A at the start to 2.82 A) while 



FIG. 4: The path of the first hydrogen in the direct diffu- 
sion path from above (top) and the side (bottom). All the 
hydrogen positions are shown, along with the initial and final 
positions of the atoms to which the hydrogen bonds (with the 
final position shown in a lighter shade). The final position of 
the ad-dimer cannot be seen in the top image since it is di- 
rectly below the initial position. Bonds (or lack of bonds) are 
produced by the imaging software, and should not be taken 
as definite indications. 



the bond from the hydrogen to the ad-dimer is more ex- 
tended (from 1.51 A at the start to 1.85 A). Inspecting 
the charge density, it is clear that the ad-dimer remains 
bonded (though weakly) and that the H has made a weak 
bond to the substrate dimer (which is 2.32 A away) as 
well as maintaining a slightly weakened bond to the ad- 
dimer. It is this lengthening and weakening of bonds 
at the saddle point that causes the high barrier. As- 
suming Arrhcnius behaviour and an attempt frequency 
of 10 13 Hz, we find a hopping rate of ~ 10 _10 scc _1 at 
450K, which is many orders of magnitude below the ob- 
served rate. 

The diffusion barrier into and out of the meta-stable 
state is also shown in Fig. |3| with open squares and a dot- 
dashed line. The barrier from the start to the meta-stable 
state is 1.24 eV, while the barrier from the meta-stable 
state to the end is 0.56 eV (and the reverse path, from 
the meta-stable state to the start is 0.66 eV). 

The pathway from the starting position to the meta- 
stable state (shown in Fig. |^) involves considerable re- 
arrangement: first, the hydrogen inserts into the bond 
between the ad-dimer (labelled 'C in Fig. |J) and the 
substrate dimer (labelled 'A' in Fig. |J); at the saddle 
point, the hydrogen is 1.65 A from the ad-dimer, and 
1.94 A from the substrate dimer (compared to an equilib- 



FIG. 5: The path of the first hydrogen from the ad-dimer to 
the meta-stable state shown in views from above (top) and the 
side (bottom). All the hydrogen positions are shown, along 
with the initial and final positions of the atoms to which the 
hydrogen bonds (with the final position shown in a lighter 
shade). The side view is shown rotated by 180° relative to 
the view from above as the image is clearer. The final position 
of the ad-dimer cannot be seen in the top image since it is di- 
rectly below the initial position. Bonds (or lack of bonds) are 
produced by the imaging software, and should not be taken 
as definite indications. 

rium distance of 1.51 A), while the distance between the 
ad-dimer and the substrate dimer is 3.10 A (compared 
to 2.48 A at the start); second, after the hydrogen has 
transferred to the substrate dimer, the ad-dimer bonds 
to a second layer atom in the substrate (labelled 'B' in 
Fig. ||) ; third, the substrate dimer bonds back to the ad- 
dimer and breaks its bond to the second layer atom in 
the substrate. The first part of this rearrangement is the 
area where most of the energy change happens: the en- 
ergy actually falls by about 0.6 eV during the second and 
third parts of the rearrangement. 

The pathway from the meta-stable state to the end po- 
sition (shown in Fig. ^) is much simpler, involving only 
the movement of the hydrogen from one end of the sub- 
strate dimer to the other, while the substrate dimer 'C 
reforms its bond to the second layer atom 'B'. At the 
saddle point, the hydrogen is 1.72 A from the substrate 
atom and 2.05 A from the end atom. 

The barrier of 1.24 eV from the starting point to the 
meta-stable state fits extremely well with the observed 
temperature behaviour: at 450K with an attempt fre- 
quency of 10 13 sec^ 1 , it would correspond to a hopping 
rate of 0.044 Hz, or one hop every 23 seconds. But this is 
only into the meta-stable state, and there are two low en- 
ergy paths out of that. The hopping rate from the meta- 
stable state to the end state is ~ 4 x 10 6 Hz, while from 
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FIG. 6: The path of the first hydrogen from the meta-stable 
state to the substrate dimer shown in views from above (top) 
and the side (bottom). All the hydrogen positions are shown, 
along with the initial and final positions of the atoms to which 
the hydrogen bonds (with the final position shown in a lighter 
shade). The side view is shown rotated by 180° relative to 
the view from above as the image is clearer. Bonds (or lack 
of bonds) are produced by the imaging software, and should 
not be taken as definite indications. 

the meta-stable state to the start state is ~ 2 x 10 5 Hz, so 
that only 10% of meta-stable states would return to the 
starting point. We also expect that the equilibrium pop- 
ulations of the start and end states would differ, since 
the end state is 0.53 eV lower than the start (roughly, 
at 450K, we would expect a relative population about 
10 6 times higher in the lower state). There is also the 
question of whether the hydrogen could return, via the 
meta-stable state, from the end to the start. The barrier 
from the end point back to the meta-stable state it is 
1.66 eV, making it extremely unlikely that the hydrogen 
would return to the meta-stable state (and even if it did, 
it would be ten times more likely to drop back to the 
end state than to return to the start state) . Clearly, it is 
the low barrier from the starting state to the meta-stable 
state that allows the first part of the dehydrogenation 
of the ad-dimer to proceed, and the energy difference be- 
tween the start and end points, as well as the high barrier 
out of the end state that makes the reaction effectively 
irreversible. 



B. The Second Hydrogen 

Once the first hydrogen has diffused off the ad-dimer, 
we retain it on the substrate, as shown in Figure |l|(b). 
This is computationally convenient (as it maintains a 



FIG. 7: A graph showing energy barriers for the diffusion of 
the second hydrogen off the ad-dimer (with the first hydrogen 
already on the substrate). The energies are given relative to 
the starting point, while the z-axis gives the distance from 
the starting point of the hydrogen. The open circles show 
direct diffusion (proceeding without the meta-stable interme- 
diate state). The open squares show the diffusion via the 
meta-stable state (whose position is marked with a vertical 
dotted line at 3. 41 A, and whose structure is shown in Fig. |^). 
The lines (long dashes for direct diffusion and dash-dotted for 
diffusion via the meta-stable state) are spline fits to the data, 
and are given as guides to the eye. 

filled set of bands) but also physically reasonable: hy- 
drogen does not begin diffusing along the dimer rows on 
Si(001) at an appreciable rate until about 550K with a 
barrier of 1.68 eVEl 

As with the first hydrogen, the second hydrogen can 
diffuse either directly, or via a meta-stable state, which 
is exactly equivalent to the meta-stable state for the first 
hydrogen (shown in Fig. |2|) and therefore not illustrated 
here. As before, we will discuss these results separately, 
starting with the direct diffusion. 

The diffusion barrier for direct diffusion is shown in 
Figure [l], plotted with open circles and dashes. The shape 
is much broader than for the first hydrogen's direct path, 
with a lower barrier of 1.59 eV. The reason for this can 
be seen in part in the atomic positions, which are shown 
in Figure |^. This is a little more confusing than previous 
plots, as the positions both of the hydrogen and the ad- 
dimer atom to which it is bonded have been plotted. As 
the hydrogen moves across towards the substrate dimer 
(the end point), the bond between the silicon atoms in 
the ad-dimer breaks, with the atom that the hydrogen 
is bonded to following the hydrogen as it diffuses. At 
the saddle point, the hydrogen is 1.67 A from the ad- 
dimer and 2.03 A from the substrate dimer, while the 
distance between ad-dimer atoms is 4.48 A. Beyond this 
point, the hydrogen transfers to the substrate, and the 
ad-dimer reforms slowly. 

The bond in the clean or partially clean ad-dimer is 
not as strong as the other bonds to the substrate, which 
explains why the energy cost for breaking it is rela- 
tively small, and why this pathway is followed in con- 
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are now identical, meaning that 50% of meta-stable states 
will return to the starting state. However, the end state 
is 0.48 eV lower in energy than the start, so that (as 
before with the first hydrogen) we would expect the pop- 
ulation in thermal equilibrium at 450K to be about 10 6 
times higher in the end state than the start state. The 
barrier from the end state back to the meta-stable state 
is 1.63 eV, which again makes the reaction effectively ir- 
reversible. Of course, the clean ad-dimer can also diffuse 
away along the. trench between dimer rows, with a bar- 
rier of 1.15 eVu, which would make the reforming of the 
hydrogenated ad-dimer impossible. 



CONCLUSIONS 



FIG. 8: The path of the second hydrogen in the direct dif- 
fusion path shown in views from above (top) and the side 
(bottom). All the hydrogen positions are shown, as are the 
positions of the atoms to which the hydrogen bonds. Bonds 
(or lack of bonds) are produced by the imaging software, and 
should not be taken as definite indications. 

trast to the pathway for the first hydrogen. Even with 
the reduced barrier, the hopping rate at 450K will be 
~4x 10~ 6 sec -1 , which is still far too low to be consis- 
tent with the experimental observations. 

The diffusion barrier into and out of the meta-stable 
state is also shown in Figure 0, plotted with open squares 
and dot-dashed lines. The barrier from the start to the 
meta-stable state is 1.14 eV, while the barrier from the 
meta-stable state to the end is 0.58 eV (and the re- 
verse path, from the meta-stable state to the start is 
also 0.58 eV). The atomic positions are almost identical 
to those for the first diffusion (shown in Figs. || and |^) 
and are not shown (though these figures, and many other 
pieces of supplementary information such as animations 
of the processes can be found elsewhereEJ) . 

The barrier of 1.14 eV is 0.1 eV lower than the bar- 
rier for the first hydrogen, suggesting that once the first 
hydrogen has diffused off the ad-dimer, the second will 
follow slightly more quickly; it is still in excellent agree- 
ment with observed experimental behaviour. The barri- 
ers from the meta-stable state to the start and end states 



We have presented ab initio calculations, modelling the 
diffusion of hydrogen off a hydrogenated ad-dimer, which 
is a key stage in gas-source growth of Si(001). We have 
shown that the diffusion proceeds via a meta-stable inter- 
mediate, and that the energy barriers calculated (1.24 eV 
for the first hydrogen and 1.14 eV for the second hy- 
drogen) are in excellent agreement with temperatures at 
which these features are observed in experiment. 

We have used the climbing image nudged elastic band 
method to find the diffusion barriers, and have found it 
to be extremely effective, particularly for the direct diffu- 
sion which was difficult to model simply by picking a sin- 
gle constraint. However, the problem of exploring phase 
space is still a difficult one, as the existence of the meta- 
stable state (which was discovered through application of 
a single constraint) shows. There are techniques fax ex- 
ploring energy surfaces, such-as the dimer methoccj and 
variants on hyperdynamicsc3, but there is still a large 
amount of work to be done in this field. 
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